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Abstract — Weak noise smooths out fractals in a chaotic 
state space and introduces a maximum attainable resolu- 
tion to its structure. The balance of noise and deterministic 
stretching/contraction in each neighborhood introduces lo- 
cal invariants of the dynamics that can be used to partition 
the state space. We study the local discrete-time evolution 
of a density in a two-dimensional hyperbolic state space, 
and use the asymptotic eigenfunctions for the noisy dynam- 
ics to formulate a new state space partition algorithm. 



1. Motivation and outline 



Chaotic systems' main feature is their high sensitivity 
to initial conditions. That makes direct numerical integra- 
tion of the equations difficult and often calls for alterna- 
tive methods for the evaluation of long-time averages of 
observables, such as decay of correlations, diffusion coef- 
ficients, energy spectra, or escape rates [ ]. To properly 
weigh these averages, one needs to understand which re- 
gions of the state space are more or less relevant for the 
dynamics, in other words make a partition [3]. Invariants 
of the dynamics such as unstable periodic oribts have been 
successfully used to partition the state space [1]. 

However, noise, modelled by stochastic variables, erases 
periodic orbits. One has to look for new invariants. For 
that reason, we previously have studied [4, 2] the evolu- 
tion of densities of trajectories and determined eigenfunc- 
tions of the local Fokker-Planck operator in the vicinity of 
the deterministic periodic orbits. The eigenfunctions are 
then used to partition the state space. All that was done in 
discrete time in one dimension. In order to develop a simi- 
lar algorithm in higher dimensions, the first step is again to 
study the evolution of densities in the neighborhood of the 
periodic points of the deterministic system. 

In the present contribution we focus on the asymptotic 
evolution in two dimensions, forward and backward in time 
of a noiseless hyperbolic map (sect. 2), to which we succes- 
sively add weak, uncorrelated, isotropic noise (sect. 3). In 
both cases the densities asymptotically align with the unsta- 
ble (stable) direction of the monodromy matrix when iter- 
ated forward (backward) in time. We finally use our results 
to propose a definition of a neighborhood for an optimal 
partition of the state space. 



2. Deterministic evolution 

We start by reviewing the deterministic evolution of den- 
sities and observables in the neighborhood of a fixed point 
x of the two-dimensional map x' = f(x). We assume that 
the fixed point is hyperbolic, i.e., that the Jacobian matrix 
evaluated at the fixed point, 



Mij(x ) = 



dfi(x) 



(i) 



has eigenvalues |A S | < 1, |AJ > 1. 
Consider the simplest example, a map 

f(x) = (A s x,A u y) 



(2) 



that is contracting along the x-axis and expanding along 
the v-axis. Now consider a density of trajectories p(x), for 
instance a Gaussian placed around the fixed point of f(x), 
and apply the Perron-Frobenius operator [1] to it: 



£p(x) = dz£(x-f(z))p(z) 



1 
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so that, after n iterations, 



£T P (x) 



1 



\A u f \A S 



(3) 



(4) 



One can see this as a density, which is losing mass by a 
factor of |AJ _1 at each iteration. This expression can be 
renormalized by a factor of \A u \ n , when taking the limit 
n — > oo. If the initial density is a normalized Gaussian 
p(x) oc exp [-(x 2 + y 2 )/2<x 2 ], we obtain 



lim \A u \ n £ 



lim 

n— >oo 



4 



(5) 
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meaning the limiting density is supported on the y-axis, the 
unstable manifold of the fixed point of the map. 

Of course this was the simplest possible example, given 
that the contracting and expanding directions of the fixed 
point are already separated by the coordinates. This is not 
the case in general, and one needs to do something differ- 
ent from what we just described. We will follow Rugh's 
formalism [8] for a general two-dimensional map f(x) with 
a hyperbolic fixed point xq: the equation f y (xi,yi) = V/ has 



a unique solution, which we can call <f> s (Xi,yf), which is an- 
alytic and a contraction. On the other hand, one can define 

M*h y/) = fx (xu (f> s (xi, y f )) (6) 

and then rewrite f , 

f (x h (f> s (x u y f )) = (<f>u(xi, y/X y/) , (7) 

in terms of the pinning coordinates (x;,y/), that is the con- 
tracting coordinate of the initial point, xi and the expanding 
coordinate of the final point, yf. It is important to remark 
that both <p u (xi,yf) and (f> s (xi,yf) are contractions on their 
supports [ ] . In particular, for fixed x t , 



lim (p n s (x h y f ) = W s (xi) 



(8) 



with W s (xi) such that (jc,-, W s (xi)) parametrizes the stable 
manifold of the map f(x). Similarly, 



lim <ft(x h y f ) = W u (y f ) 



(9) 



where (w u (yf),yA defines the unstable manifold of f(x). 

Now we will obtain again the limit (5) for the evolution 
of a density carried by the Perron-Frobenius operator. We 
will study the evolution (4) inside a space average of an 
observable a(x). 

< a >n = TT7, f dxa(x)[£ n p(x)] (10) 
|VH| J M 

As we will show, the support of the observable inside the 
average corresponds to the support of the mapped density. 
Eq. (10) is equivalent to a more familiar expression for the 
space average, which we can easily write by letting X" act 
on its left on the observable a, in which case it becomes the 
Koopman operator 7C [1], 

f dxa(x)[£ n p(x)] = f dxYK n a(x)]p(x) 

JM JM 



f. 

JM 



dxa(f n (x))p(x) (11) 



We now change the coordinates in the last integral accord- 
ing to the transformation (7): 

J dx t dy f a (f u (x h y f ), y f ) p (x u f s (*,-, y f j) det (d 2 </> n s (x h y f )) 

(12) 

where / is the new domain of integration and 
det (d2<p"(xi,yf)) is the determinant of the Jacobian 
of the change of coordinates in the integral. For n — > oo 
one gets, according to (8) and (9), 

j dx i dy f a(w u (yf),y f )p(x i ,W s (x i ))dtt (d 2 ^(x h y f )) . 

(13) 

We can see that the observable a ends up being supported 
on the unstable manifold of the map. Intuitively, the ini- 
tial observable stretches and contracts respectively along 



the unstable and stable manifolds, so that it asymptotically 
survives on the only region of the state space (the unstable 
manifold) where it cannot be crushed by the contraction. In 
the separable case (2), the average (10) is proportional to 

J dxdy a(x, y)£ n p(x, y) — > j^-p- J dxdya(x, y)S(x) = 

j^r fdya(0,y). (14) 

In other words, knowing the support of the observable a 
inside the average is equivalent to knowing the support of 
the time-forward evolved density p (in this case the y-axis, 
cf. (5)). 

Now we consider the time-backward evolution, de- 
scribed by the Koopman operator 7C: 

f M dxa(x)[K n p(x)]^ 
ldxidy f a{x u W s (xi))p(w u (y f \ y f )dct (d 2 ^(x h y f p) 

where we switched to pinning coordinates as in (13), with 
p and a inverted with respect to the previous case, and then 
took the limit n — > oo. This time the observable a is asymp- 
totically supported on the stable manifold. 

Analogous results are found for a periodic point x a of a 
map f(x), for it can be regarded as the fixed point of the 
iterated map f^(x), with n p period of the cycle. The time- 
forward (backward) evolution aligns observables and thus 
densities to the unstable (stable) eigenvector of the mon- 
odromy matrix 



dx; 



(16) 



evaluated at the periodic point x a 



3. Adding noise 



We now add weak noise to the map f(x). In the vicinity 
of any point x a , Gaussian densities are mapped forward in 
time by the Fokker-Planck operator £ Fp [2, 6] 

r. (ry \ — C U 7 ^ 0-\{Za+\-M a Z a ) T T{Za+\-M a Z a )-Za T TT Za 

Pa+lKZa+l) ~ I YaZa\e 2 ® a 

1 1 _ T 1 _ 

(17) 



C, 



Qa + \ • 



a+\ 



where we defined local coordinates z a = x — x a . Here the 
noise is described by the symmetric and positive definite 
diffusion tensor A. If the density is a Gaussian distribu- 
tion, we can recast the problem in terms of its covariance 
matrix [2]: 

Q a+ i =M a Q a M a T +A. (18) 

The long-time limit is given by the fixed-point condition 
Qa = Qa+i, valid when the dynamics is contracting (M 
has all eigenvalues |A/| < 1). This condition states that the 
covariance matrix must be invariant under the combined 
action of the deterministic contraction and expansion by 



weak noise after one time step. Let S be the matrix which 
diagonalizes M. If we make the further transformations 
Q -> S- 1 Q(S~ 1 ) T = Q and A -> 5" 1 A(s- 1 ) t ee A, the 
solution to the fixed-point condition for contracting maps 
reads [2] 

1 



1-A,A,- 



(19) 



We obtain the time-backward evolution by taking the ad- 
joint of the operator in (17). Like before, an equation is 
derived for the mapping of the co variance matrix: 



M a Q a M a T = Q a+1 + A. 



(20) 



If the deterministic dynamics is expanding (M has all 
eigenvalues |A;| > 1), (20) becomes a fixed point condition 
by setting Q a = Q a +\. We find that 



Qu = 



(21) 



where we applied the same diagonalization transformation 
by means of the matrix S . 

Typically in a chaotic system the matrix M has both con- 
tracting and expanding directions, so that neither the solu- 
tion given by (19) nor by (21) applies. In what follows we 
study the evolution of the covariant matrix both forward 
and backward in time, looking for an asymptotic limit. We 
start forward in time, iterating (18) in the neighborhood of 
a fixed point of f(x): 

Q n = A + MAM T + M 2 A(M 2 ) T + • • • + M n Q (M n ) T . (22) 

As in the deterministic case, let us first see how the g's 
map when M and the initial go are diagonal: each matrix 
element Qu obeys the sum (22), which diverges for |A/| > 1, 
so that Q^ 1 vanishes. On the other hand, it converges to 
Qu — > Au/(l - A 2 ) [4] if |A;| < 1. As a result, the axes of 
the Gaussian e~ zTQ z asymptotically survive in the stable 
directions only, while the whole density is supported along 
the unstable directions, like in the deterministic separable 
case (cf. (5)). In two dimensions, the asymptotic density 
looks like the Gaussian-shaped tube in figure 1(d). 

We will now explain the asymptotic evolution of the axes 
of the ellipsoid the Gaussian is supported on when M is not 
diagonal, mimicking an argument of Ott's [5]. First of all, 
if the noise is isotropic, A(x) - 2D1, the diffusion tensor 
A is replaced by a scalar diffusion constant D that can be 
factored from the sum (22). Consider a vector v: 



(a) 




(c) 



Figure 1 : Evolution of a Gaussian density centered at the 
fixed point of a hyperbolic map: (a) initial bump; (b) after 
one iteration; (c) after two iterations; (d) asymptotic limit. 



us now write M n = SA n S l , with S = (e\, e U9 e^), 
and S~ l = (e[, e' u , e' N ) T . So now 



SA"(e' u -e[,...,l,...,e' u -e' N y = 
S{K n x e[-e' u ,...,K,...,K n N e' N -e' u ) T 



(24) 



When n — > oo the uth component of the vector weighs 
above all others, making the result proportional to 



S(0,...,1,...,0) T = e u . 



(25) 



Thus any vector is eventually stretched and rotated toward 
the most unstable direction of M. It is straightforward 
to show that the last term in the sum (22), M n Q (M n ) T , 
asymptotically behaves likewise. Consequently, goo and 
thus g" 1 are also aligned with the most unstable eigen- 
vector of the monodromy matrix M. Numerics (figure 1) 
help us visualize the result in two dimensions, where M has 
one stable and one unstable directions: an initial isotropic 
Gaussian develops into a 'tube', infinitely extended along 
the unstable manifold of M, and having a Gaussian section 
in the orthogonal direction, due to the balance of noise and 
deterministic contraction. 

One can repeat the above reasoning when applying the 
adjoint Fokker-Planck operator £} Fp . In this case we in- 
vert (20) (the unknown being Q a ), and iterate it n times to 
get 



M n (M n ) T v = M n Y, ajA n /j - M n a u A n u e' u (23) Q-n = M~ n Q (M- n ) T + M~ n A(M~ n ) T + ... + M" 1 A(M" 1 ) T . 



where we first wrote v in terms of the eigenvectors e'. of 
M T , then we applied (M n ) T to the same vector and finally 
we observed that the result asymptotically aligns to the di- 
rection e' u of the most unstable eigenvalue of M T , A u . Let 



(26) 

This is similar to (22), except the monodromy matrix is in- 
verted, so that stable and unstable eigenvalues are swapped, 
and the argument (23)-(25) results in any vector being 
stretched and rotated toward the most stable direction of 
M. 



The observations made for a fixed point of the map can 
be extended to a periodic orbit of arbitrary period n p . We 
start again from (22): 

Q a = A + M._ 1 A(M«_ 1 ) T + • • • + M\Q a . np (M n /_ np ) T . (27) 
We define 

Ap^ = A + M«_ 1 A(M a _ 1 ) T + 

and then the asymptotic evolution of the covariance matrix 
Q around the periodic point x a can be written as 

Q P ,a = M p , a Q a M p , a T + A p , a , (29) 

where M PA - M n a p , the latter defined in (16). The problem 
reduces to the previous case, which shows that the leading 
eigenvector of the asymptotic covariance matrix (or the ma- 
jor axis of the ellipsoid) aligns to the most unstable (stable) 
eigenvector of M Pt0 in the time-forward (backward) evolu- 
tion. 

Following the technique we used in one dimension [ ], 
we define the neighborhood M a of the periodic point x a as 
the intersection of the supports (within a Icr confidence) of 
the ground-state local eigenfunctions of £ FP and £ Fp . We 
use these regions to cover the non-wandering set of the sys- 
tem, starting with the periodic orbits of the shortest period, 
and increasing the period until neighborhoods significantly 
overlap (see figure 2(d)). 

4. Summary 

We have studied the asymptotic evolution of Gaussian 
densities of trajectories in the neighborhoods of hyperbolic 
periodic points, first in a deterministic map and then in the 
presence of weak, uncorrected, isotropic noise. We inves- 
tigated both time-forward and -backward dynamics. The 
latter was realized by means of the adjoint of the evolu- 
tion operator. As it turns out, the densities asymptotically 
align with the most unstable (stable) direction of the mon- 
odromy matrix when iterated forward (backward) in time. 
Using both asymptotic densities, we proposed a definition 
for a neighborhood which should then be used to partition 
the state space. 
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